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1. Introduction 

The 3D Ising model has, over the last 25 years, been used to test the accuracy of various analytical and 
numerical methods for solving Statistical Mechanics systems. In 1992 we presented results of simulations 
on 64'^ and 128'^ lattices using the Monte Carlo Renormalization Group (MCRG) method [Q. While that 
work improved on previous MCRG estimates 1^ [||), it left us with four unanswered questions. The first 

and most tantalizing was are the exponents rational numbers, i.e. v = 0.625 and t] = 0.025. Second, 

a more precise determination of the corrections-to-scaling exponent is needed as it is the largest source of 
systematic errors. Third, we wanted to resolve/understand the differences between our MCRG and existing 
finite size scaling/ e-expansion results. Lastly, we wanted to investigate whether hyperscaling holds for this 
model. This talk is a summary of the current status of our calculations. 

In order to address these issues we have extended the calculations in the following ways. We have made 
higher statistics runs on 64^, 128^, and 256^ lattices at K = 0.221652 and 0.221655. On these lattices 
we have evaluated, in addition to correlation functions needed for MCRG studies, quantities needed for 
finite size scaling (FSS) analysis and the calculation of the renormalized coupling gn. As a result we have 
better estimates of the critical coupling the exponents v and 77 from both MCRG and finite-size scaling 
(including histogram re-weighting) analysis, and can address the issue of hyperscaling violations. Our new 
estimate of the corrections-to-scaling exponent lo « 0.7 is significantly smaller than that from other methods. 

All the simulations have been done on the Thinking Machines CM-5 (at the ACL at LANL) and CM-5E 
(at TMC) computers. We used the Swendsen-Wang cluster update algorithm [|| and a 250-long 64-bit wide 
shift-register (Kirkpatrick-StoU) random number generator in each vector unit. The new results agree with 
our previous calculation and those in |^ and |^ . Each of these calculations used a different random number 
generator, so their consistency suggests that there is no obvious bias in the sequence of random numbers 
generated (see P. Coddington's talk on random number generators at this workshop). Our most extensive 
results are at Knn — 0.221655, which is our present best estimate of Kc, and the statistical sample consists 
of 600K, 500K, and 400K measurements on 64^, 128^, and 256^ lattices respectively. 

The details of our implementation of the MCRG method are the same as in [|l|. The only change is 
that we have added 3 more even couplings (for a total of 56) and one more odd couplings (total 47). The 
original 53 even and 46 odd couplings were contained in either a 3 x 3 square or a 2^ template [§. The new 
couplings are those obtained by adding a fourth spin along the cartesian axis to the 3x3 template. 

We store the magnetization and energy for each configuration, from which we can calculate quantities 
like the specific heat, susceptibility. Binder's cumulant U — 3 — {ra^) / {m?)"^ , etc.. These results are then 
evaluated as a function of K, in a small neighborhood of K simulation, using the histogram re- weighting 
method ||]. The finite size analysis of these quantities follows the work of Ferrenberg and Landau Q, i.e. 
without corrections to scaling terms. To calculate gn, we also need the finite lattice correlation length ^. 
This is calculated in two ways: 

(1.1) 

k-^\s{kY ) ^ ' 

where Sik) = J2x y z ^i^^V^ ^)^^'''^ ■ investigate the 5 lowest momenta but present results only for the 
lowest, kz = 2tt/L, as it has the best signal. With ^ in hand we calculate ga defined as |Q 

\ (m2)2 
This is expected to scale as 

gniK.L) ~ * (1.3) 

where w* = + dv — 2A)/j/ is referred to as the anomalous dimension of the vacuum. If hyperscaling holds 
then g'^ — > finite non-zero constant as L ^ 00 and K —>■ K*. 

For the purpose of error analysis the data have been divided into bins of size 10,000 measurements. All 
errors are then calculated by a single elimination jackknife procedure over these bins. This talk is organized 
as follows. We first summarize the MCRG results, then compare them with FSS estimates, and finally give 
data for gu, both at Ksim — 0.221655 and along the MCRG flow, and discuss hyperscaling. 



gn{K,L) = ( 3 - £^ ) . (1.2) 
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2. Nearest-Neighbor Critical Coupling K^^ 

We calculate i^^^ using the two lattice MCRG method. The results for the two starting couplings 
Knn — 0.221652 and 0.221655 and for the two pairs of lattices, 256^ versus 128'^ and 128'^ versus 64^, are 
shown in Table where for ease of comparison we have also reproduced results from We expect the 
deviations K'^{oo) — K'^{n), where n is the blocking step, to converge by the geometric factor A„/At as a 
function of the blocking level n. For the 3D Ising model Xu/Xt ~ 6, and our data for (see Table |l|) do 
roughly show convergence by this factor. 

For both starting couplings the data show that the new estimates from the 128'^/64'^ lattices analysis 
are consistent with our earlier results and give = 0.221652. However, the 256^/128^ resuhs show better 
convergence with respect to blocking steps, and also show a systematic shift at all blocking levels compared 
to the 128^/643 data. Our new estimate from 256^/128^ comparison is = 0.221655 ± 0.000001 where 
only the statistical error has been quoted. This shift in with lattice size shows that there could still be 
finite size corrections at the level of statistical errors, i.e. a systematic error of +0.000001. So our final best 
estimate is 

= 0.221655 ± 0.000001 +^ °"™°' , (2.1) 

which is consistent with the recent result 0.2216546(10) obtained using FSS analysis by Blote et al. In 
view of this we take the results at fC„„ = 0.221655 to represent the critical point values. 



Level 


0.2216550 


0.2216520 


0.2216550 


0.2216540 


0.2216520 


0.2216440 


n/m 


256^/1283 


256^/1283 


128^643 


128^/643 


1283/64=' 


128^/643 




New 


New 


New 


§ 


New 


i 


2/1 


217± 11 


184 ± 15 


095 ± 13 


070 ± 16 


046 ± 23 


095 ± 17 


3/2 


469 ±9 


458 ± 18 


413 ± 15 


406 ± 18 


394 ±27 


417 ±22 


4/3 


537 ± 10 


534 ± 21 


504 ± 16 


500 ±21 


501 ±31 


508 ± 24 


5/4 


549 ± 10 


549 ± 22 


523 ± 16 


514 ±26 


516 ±37 




6/5 


547 ± 10 













Table 1. Estimates of -f^^„ as a function of the blocking level and K simulation- For brevity only the last 
three decimal places have been quoted, so 547 ± 10 is short for 0.2216547 ± 0.0000010. We have 
included our old data from to facilitate comparison. The highest blocked lattice is 4^ except at 
Knn = 0.221652 where it is only 8^ on 128^ and 256^ lattices. The quoted errors are the statistical 
errors after averaging the data over the 56 even operators. The data show a systematic shift between 
the 2563/128^ and 128^/643 lattices. 



3. Correlation length exponent 

The correlation length exponent v is determined from the leading even eigenvalue At of the linearized 
transformation matrix 7^^, 

_ 1 In At 

yt ^ - = -—r , (3.1) 
v ino 

where b — 2 is the scale factor of the majority rule blocking transformation. Our preferred data for At 
{Ksim — 0.221655) is shown in Table || as a function of the blocking step and lattice size. There are three 
possible sources of systematic errors that affect the L oo and K — + Kc estimates for At. These are 
1. The number of operators measured, i.e. the truncation errors in evaluating eigenvalues from a finite 
dimensional 72^^. We find that the number of operators needed to achieve convergence increases with 
the blocking level n. The data show that with the 56 even operators the eigenvalues show convergence at 
all levels. (The same is true in the sector of odd interactions from which we extract the exponent rj) . Even 
at the highest blocked level there is no detectable variation after including 30 operators. Unfortunately, 
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the convergence with the number of operators is not monotonic and there is no independent way of 
confirming that the results have converged. Thus, we cannot estimate the possible error due to lack of 
convergence, and guess that it is smaller than the statistical error. 

2. Finite size effects on blocked lattices. It has been observed in jl], and Q that finite size effects are 
discernible only when blocking from 8'^ 4^ lattices or smaller. The correction increases the estimate 
of At. In Ref. the correction in Xf^^ was estimated to be 0.02. Our estimate based on comparing 
256'^, 128^, 64^ lattices is « 0.01. Applying +0.01 as the correction to our 256^ lattices data, we get the 
lower limit A^^'* = 3.008 corresponding to ly ~ 0.6294. We discuss the L ^ oo limit below. 

3. Error in the estimate of K^„. The value of Af also increases with K^^ as shown by the data in Fig. |^, 
and in Ref. . The dependence of Ai on K^^ is marginal on the first couple of blocking levels and then 
increases rapidly with n for n > 3. Since our estimate of -R'^„ is converging from below, our results at 
K = 0.221655 may underestimate A*. 

The bottom line is that the systematic effects discussed in items 2 and 3 will tend to increase At or equivalently 
decrease v. 





At 


A,, 




256^ 


128=^ 


64^ 


256^ 


128^ 


64^ 


0/1 


2.681(2) 


2.685(2) 


2.684(3) 


5.4948(06) 


5.4941(06) 


5.4948(08) 


1/2 


2.843(2) 


2.847(2) 


2.847(2) 


5.5050(02) 


5.5052(02) 


5.5050(03) 


2/3 


2.930(3) 


2.930(3) 


2.930(3) 


5.5501(02) 


5.5499(03) 


5.5494(07) 


3/4 


2.969(3) 


2.973(4) 


2.971(5) 


5.5741(04) 


5.5745(08) 


5.5701(21) 


4/5 


2.995(5) 


2.985(7) 




5.5845(11) 


5.5826(31) 




5/6 


2.998(7) 






5.5850(31) 







Table 2. Estimates of At and Xh as a function of the blocking level and lattice size for K — 0.221655. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Fig. 1: Estimates of At as a function of the blocking level n for simulations at different 
couplings K . The best estimate of the extrapolated value with ui = 0.70 (uj = 0.85 ) 
is from the linear fit to the diamonds skipping n = 0, 1 (n = 0) points. The fit to 
all the diamond points is, for comparison, shown by the dotted line. 
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Finally, we are interested in the value of At at the fixed point. To obtain this we extrapolate At versus 
the blocking level n using 

At(n) = A* + atfe-"" . (3.2) 

where lu — 9/i> is the leading correction-to-scaling exponent. There are two issues that need to be addressed 
in doing this extrapolation in the number of blocking steps n (i.e. the L — > cxd limit). The first is the value 
of u! and the second is whether the fit should exclude the first few blocking steps to avoid transients, to 
account for which requires further corrections to the leading behavior shown in Eq.( ^.2| ). The calculation of 
LU is discussed in the next section and our present estimate uj = 0.7 is surprisingly low. We, therefore present 
an analysis for uj — 0.7 and 0.85, where the second estimate is roughly what is given by other methods (FSS, 
e-expansion, etc. See ^ for a very recent survey). The question of transients is completely empirical, i.e. 
we neglect data at initial blocking steps until 1. 

On basis of the quality of the fit to the K = 0.221655, L ~ 256 data the best estimates for the two 
extreme values of uj are 

A* = 3.028(3) =^ ly = 0.6256( 5) (u 0.85, n = 2 - 6) , 

(3.3) 

A* = 3.033(6) =^ ly ^ 0.6247(10) (to = 0.70, n = 3 - 6) . 

These two estimates are consistent, we therefore take the mean value and the larger of the two errors to get 
our present best estimate ly = 0.625(1). To improve this result will require a better estimate of uj and data 
on larger lattices (more blocking steps). 



4. Correlation function exponent rj 

The correlation function exponent 77 is given by 

rj = d + = d + 2-2vh , (4.1) 

logfe 

where 6 = 2, d = 3 and A/i is the largest eigenvalue of 7^/3 constructed from the odd interactions. The 
discussion of the type and sign of the various systematic errors in the extraction of \h is the same as for At . 
The raw data are shown in Table ^, and the value of finite size correction we apply to A^^'' is 0.002. Then, 
from the L = 256 data (A^^'* = 5.587(3)) we get the upper bound rj = 0.0359(16). 

To extrapolate to L — > cxd we proceed in exactly the same way as for At. However, as exemplified by 
Fig. H, the points at n = 4, 5 show significant deviations from the linear fits. Even though the fit with 
UJ — 0.85 is somewhat better, the n — 5 point raises questions about the validity of the linear extrapolation. 
There are two possibilities. One, the value flattens out at Xh = 5.59, in which case 77 — 0.034. Second, the 
points at higher blocking levels are not well determined (note the large dependence of A^ on K at higher 
n as shown in Fig. ||. The systematic errors could therefore be larger than the statistical), and the linear 
extrapolation is valid. In the latter case one gets 

A^ 5.603(4) =^ ?/ = 0.028(2) (uj = 0.85, n = 2 - 5) , 

(4.2) 

a;; = 5.610(5) =^ T] = 0.024(3) (uj = 0.70, n = 2 - 5) . 

For errors we have used the difference between the extrapolated values with fits to n = 1 — 5 and n = 2 — 5 
points. The bottom line is that even though we have improved the estimates for A/i on lattices of size up to 
L = 256, there is still a large ambiguity in the determination of rj due to the extrapolation to L = 00. 



5. Corrections-to-scaling exponent uj 

It should be clear from the above discussion that a precise estimate of uj is very important in order to 
take the L — > 00 limit. In a MCRG calculation the correction-to-scaling exponent is determined from the 
sub- leading eigenvalue Af^2 the even sector; uj = —yt,2 — ~ log At, 2/ log 6. Only if uj is known can the 
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Fig. 2: (A) Estimates of Xh as a function of the blocking level n for simulations on 256'^ 
lattices at K ~ 0.221655. The same data is plotted for uj = 0.70 and 0.85, and 
for each case we show a linear fit to n = 2 — 5 points. The n = 5 point deviates 
significantly from the fits. (B) Xh as a function of K and n. 



exponents v and ?/, calculated along a critical RG flow, be extrapolated to the fixed point using Eq. ( [3.2| ). 
We have therefore spent considerable effort in estimating With our best data, shown in Table ^, we have 
overcome the statistical problem that plagued the data in Q. The second and third eigenvalue no longer 
merge into a complex pair when the number of operators is > 15. However, we now find that the value of uj 
decreases both with the number of operators and blocking steps; it starts off at « 0.85 on the first blocking 
step and with 10 operators, but finally settles down to 

Lu - 0.70 =^ e = ujvr^Q.U . (5.1) 

This value is significantly smaller than the world average, 6 — 0.54(3), of estimates obtained using other 
methods as presented in pi . Clearly, this issue requires further attention. 





10 Ops 


20 Ops 


30 Ops 


40 Ops 


50 Ops 


1 


0.90(3) 


1.01(4) 


0.97(4) 


0.98(4) 


0.97(4) 


2 


0.85(3) 


0.83(4) 


0.79(4) 


0.80(4) 


0.79(4) 


3 


0.81(3) 


0.76(4) 


0.73(4) 


0.72(4) 


0.73(4) 


4 


0.79(3) 


0.74(4) 


0.70(4) 


0.67(4) 


0.67(4) 


5 


0.75(3) 


0.73(3) 


0.70(3) 


0.70(3) 


0.70(3) 



Table 3. Estimates of lo as a function of the blocking level and the number of operators used in constructing 
the transformation matrix Tap. The data is for K = 0.221655. 



6. Finite size scaling using the histogram re-weighting method. 

We use the histogram re-weighting method |^ to estimate (i) the position of the maximum for various 
thermodynamic quantities, (ii) the value at this maximum, and (iii) the value of K at the point of crossing 



5 



of U and gu for two different size lattices. The method consists of building a histogram H{E,m), i.e. 
the number of configurations with energy E and magnetization m, using an equilibrium (canonical) Monte 
Carlo simulation at temperature Kgim- With this histogram, the equilibrium probability distribution at 
other temperatures K is 

H{E,m)eKY>[^KE] 



PKiE,m) = 



Y.^H{E,m)eM^KEy 



(6.1) 



where IS.K = K — Ksim- The average value of any function of E and m, Q(E,m), at coupling K is then 
given by 

{QK{E,m)) = Y,QiE,m)PK{E,m). (6.2) 

E,m 

The value of thermodynamic derivatives with respect to K are obtained from Monte Carlo measurements of 
correlation functions, 

d{Q) 



dK 



{QE)-{Q){E). 



(6.3) 



Again, by using the re- weighting technique these correlation functions can be evaluated at all temperatures in 
a certain neighborhood of Ksim- Thus, the location and magnitude of the peaks, and the points of crossings 
can be obtained from simulations at a single temperature. 

The propagation of errors under this re-weighting is not straightforward and has been dealt with by 
Ferrenberg and by Swendsen in their talks at this meeting. In our current analysis the error estimates are the 
naive statistical ones and ignore all correlations and uncertainty in determining H{E^m). We only present 
results for Kgim = 0.221655 as these have higher statistics and correspond to our estimate of the infinite 
volume Kc. With this and the re- weighted data generated from it in hand we use the appropriate finite size 
scaling relations to derive estimates of the critical exponents and temperature. 



20000 
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2000 
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T 1 \ — r 




0.6348(48) (U) 
0.6298(27) (log |m|) 
0.6295(27) (log m^) 



0.02 



0.01 



0.005 



0.002 



0.001 



y/v (x) = 1.9754(64) 



j/v (x')= 1.9719(41). 



50 



100 



200 



50 



100 



200 



Fig. 3: (A) Estimates of v from the Enite size scaling of the maxima of the derivatives of 
U, log TO and \ogm^ with respect to K. (B) The slope of log x versus logL gives 
the exponent ratio ^ jv assuming the leading order FSS relation x L'^l'^ . 



Estimate of v: This is obtained from the finite size scaling of the maxima of thermodynamic derivatives of 
J7, log|TO|, and logm^ For example, d{logm^) /dK = {jn^ E) / {m?) — {E) is calculated as a function of 
K using the histogram H{E, to) to re-weight the data. FSS analysis to extract v is done keeping only the 
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leading term in the scaling behavior 



dK 



(6.4) 



as we cannot reliably include correction terms with data at only 3 values of L. Linear fits to the maximum 
value versus L^^^ are shown in Fig. |[ The quality of the fits is exceedingly good, and the final results are 



0.6348(48) U 
0.6298(27) loglml 
0.6295(27) logm^ 



(6.5) 



The values obtained from the derivative of loglmj and logm^ agree, while that from U is higher by its Icr 
error estimate. These estimates are higher than the MCRG value by roughly one combined a. 

0.2222 



0.222 



0.2218 



0.2216 



0.2214 



"1 — 1~ 

K„= 

K^= 

K,= 0. 

K„= 0. 

-K,= 


K,= 

K = 



T 1 1 1 1 1 1 1~ 

2216508(43) 
2216543(16) |m| 
22165543(8) x' 
2216523(12) xl-log m 
2216525(15) |m|-log m' 
2216556(13) log |m| 
2216558(15) log m' 
2216554(17) U 




256^ 128" 

J I I I L 



J I I L 



J_ 



64-' 

J I L_ 



0.0005 



0.001 



0.0015 



L-v- (zy=0.625) 



Fig. 4: Estimate of Kc{L) obtained from the FSS analysis of different tliermodynamic 
quantities assuming v — 0.625. Tie data are for Ksim = 0.221655. 



Critical Coupling Kc- Kc has been calculated in two different ways. One, we considered the location 
of the maxima of different thermodynamic derivatives as a function of system size. The finite size scaling 
behavior of these Kc{L) is 

K,{L) =K, + a'L-i/"(l + + ...), (6.6) 

where the values of a' and b' are different for each thermodynamic quantity. The three data points for 
each quantity shown in Fig. ^ have been joined together by straight lines to highlight the deviations from 
linearity. With just three lattice sizes we cannot include the correction term (i.e. uj), to take into account the 
deviations from linearity. In all quantities, the magnitude of these deviations is surprisingly large considering 
that the data of Ferrenberg and Landau Q show a good fit to the linear form at smaller (L = 24 — 96) 
lattice sizes. We do not understand this discrepancy at present other than mention the possibility that in 
our 64^^ data AK may be too large for the re-weighting to be reliable. Since x' ^nd m converge from above 
while log |m| and logm^ approach it from below, we have constructed linear combinations, x' ^ log |m| and 
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|m| — logTO^, that show much smaller corrections. As pointed out in the quantity farthest away from 
is and it gives Kc with the largest errors. 

Fig. ^ was generated using v = 0.625. To determine Kc we show, in Table ^ the variation with v of 
Kc obtained using the data from 128'^ and 256'^ lattices to extrapolate to L = oo. We find that there is no 
significant variation with v in the canonical range of values 0.62 — 0.635. Different observables give Kc in 
the range 0.221653 — 0.221656 with a typical statistical error of 0.000002. For our best estimates we take 
the mean of the various estimates (excluding that from C^) with v = 0.625. The result is 

Kc ^ 0.2216544(20) (15) (6.7) 

where the first error is statistical and the second is an estimate of the systematic error based on the variation 
with observable type. This value is lower than that obtained by Ferrenberg and Landau but consistent 
with our MCRG result. 



Observable 


V = 0.62 


V = 0.625 


V = 0.63 


V = 0.635 


\m\ 


0.2216559(22) 


0.2216553(22) 


0.2216551(22) 


0.2216548(23) 


U 


0.2216538(24) 


0.2216539(24) 


0.2216541(24) 


0.2216543(24) 


log m| 


0.2216547(17) 


0.2216547(18) 


0.2216548(18) 


0.2216548(18) 


log 


0.2216545(21) 


0.2216545(21) 


0.2216546(21) 


0.2216547(21) 


Cy 


0.2216544(58) 


0.2216540(58) 


0.2216536(58) 


0.2216531(59) 


x' 


0.2216555(10) 


0.2216554(11) 


0.2216553(11) 


0.2216553(11) 


X' - log |m| 


0.2216540(17) 


0.2216540(17) 


0.2216540(18) 


0.2216540(18) 


\m\ — log m? 


0.2216529(21) 


0.2216529(21) 


0.2216529(21) 


0.2216529(21) 



Table 4. FSS estimates of K^„ from different observables as a function of v. 



A second estimate of Kc is given by the point of crossing of U{L) and gR{L) calculated on two lattices 
of different sizes. Again, the use of re-weighting trick to extend the data to temperatures in the vicinity of 
Ksim is essential. Our results are shown in Fig. g. The final values, taken from the crossing point of L = 128 
and L = 256 lattices data, are 

K^ = 0.2216560 U{Kc) 1.409(9) , 

ifS" = 0.2216551 gaiKc) = 5.23(10) . ^^'^^ 

Note that the error estimates do not take into account the correlations in the data generated by the re- 
weighting. It is interesting to note that the estimate of crossing point obtained from L = 64 and L = 128 
lattices is ~ 0.221652, indicating a convergence from below. Also, the estimates of Kc from the two pairs of 
lattice sizes are in very good agreement with the corresponding results obtained from MCRG analysis. 

Critical Exponent 7: We estimate 7 from the finite size scaling of the susceptibility. A linear fit to logx 
versus \ogL is shown in Fig. ^ for data at K = 0.221655 and for L= 64^, 128^ and 256^ lattices. From the 
slope we obtain 

7/1/ = 1.9754(64). (6.9) 

(To allow detailed comparison with the results in |^ we also give the value j/i' = 1.9719(41) obtained using 
x'.) Using the hyperscaling relation rj — 2 — ^/v we now get 77 = 0.0246(64). This estimate is consistent with 
the MCRG result obtained assuming the validity of the linear approximation. We thus quote rj = 0.025(6) 
as our best value since it covers the various estimates. 



7. Renormalized coupling gj^ and hyperscaling 

It has been pointed out by Baker and Kawashima ||l^ that the two limits L 00 and K Kc do 
not commute and that there is a discontinuity in the value of gfi when calculated in the following two 
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0.221652 0.221654 0.221656 0.221658 0.22166 0.22165 0.221655 0.22166 

Fig. 5: Plot of U and gn as a function of K obtained from the histogram re-weighting 
method. The Kc{L) estimated from the crossing point of 128 and 256 lattices data 
is indicated. Note that the values and error estimates are highly correlated. 

ways. Our approach, which is to Gx K = K'^{oo) and then slowly take the thermodynamic limit versus the 
"right" way which is to always keep L large and fixed (the precise value depends on the model and should 
be representative of the thermodynamic limit) while taking K K'^. Note that the field-theoretic value 
gR = 23.73(2) is an estimate representative of the "right" order of limits. Nevertheless, one can establish 
hyperscaling by our non-perturbative approach if we can show that converges to a non-zero lower bound 
as L — > oo. 

The data for gji in the vicinity of K'^ is shown in Fig. ^. One expects an increase in the slope with 
lattice size as the discontinuity due to the interchange of limits becomes sharper with lattice size. This is 
borne out by the data. On the other hand the crossing value shows a small decrease between the 128'^/64'^ 
and 256'^/128'^ lattices, i.e. the data does not converge from below. Based on our data (2 crossing points) 
we guess that gfi{oo) « 5. We regard this non-zero result as a weaker verification of hyperscaling than one 
would have liked. 

Another necessary condition to test whether hyperscaling holds, on basis of the data for gu in Fig. |^, is 
whether it is representative of the fixed point value. To check this we have also calculated ga on the blocked 
lattices. The data on 256^ 128^ ^ 64^ 32^ lattices is virtually identical, indicating that our estimate 
does not change along the flow to the fixed point. On smaller lattices the two methods for calculating ^ give 
different results and the use of finite lattice versus continuum energy-momentum dispersion relation makes 
a difference. We therefore consider our data on lattices smaller than 32^^ unreliable for the purposes of this 
test. 



8. Conclusions 

Our new results if^„ = 0.221655(1), v — 0.625(1), and 77 = 0.025(6) are an improvement over the 
previous MCRG values. These values continue to support the notion that the exponents are rational. We 
have also reconciled the disagreement between finite size scaling and MCRG results for K"^^ and 77 by a 
comparative study using the same data. The value of the exponent v from the two methods, however, shows 
a significant difference. One possible explanation is the lack of various corrections-to-scaling terms in our 
FSS analysis. 

The big surprise of the current MCRG analysis is the result lo w 0.7, which is significantly lower than 
all previous estimates. This issue clearly needs to be investigated further. 
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The convergence of gjj, defined to be the crossing point value in the hmit L ^ oo, seems to be from 
above. Thus, our data does not provide the desired lower bound to validate hyperscaling. We estimate 
5|j(L ~ oo) from data at the two crossing points to be ^ 5. If this non-zero value withstands further 
scrutiny, then we will have established that hyperscaling holds for the 3D Ising model. 
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